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ABSTRACT 

Aims. X-ray bursting neutron stars in low mass X-ray binaries constitute an appropriate source class to constrain masses and radii of 
neutron stars, but a sufficiently extended set of corresponding model atmospheres is necessary for these investigations. 
Methods. We computed such a set of model atmospheres and emergent spectra in a plane-parallel, hydrostatic, and LTE approxima- 
tion with Compton scattering taken into account. 

Results. The models were calculated for six different chemical compositions: pure hydrogen and pure helium atmospheres, and at- 
mospheres with solar mix of hydrogen and helium, and various heavy element abundances Z = 1, 0.3, 0.1, and 0.01 Z Q . For each 
chemical composition the models are computed for three values of surface gravity, logg =14.0, 14.3, and 14.6, and for 20 values of 
the luminosity in units of the Eddington luminosity, L/Ledd, in the range 0.001-0.98. The emergent spectra of all models are redshifted 
and fitted by a diluted blackbody in the RXTE/PCA 3-20 keV energy band, and corresponding values of the color correction (hardness 
factors) f c are presented. 

Conclusions. Theoretical dependences / C -L/Ledd can fitted to the observed dependence K~ >/4 -F of the blackbody normalization K 
on flux during cooling stages of X-ray bursts to determine the Eddington flux and the ratio of the apparent neutron star radius to the 
source distance. If the distance is known, these parameters can be transformed to the constraints on neutron star mass and radius. The 
theoretical atmosphere spectra can also be used for direct comparison with the observed X-ray burst spectra. 
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1. Introduction this case the photosphere radius rapidly increases and the ef- 

fective temperature decreases allowing th eir easy identification 
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Neutron stars (NSs) are the densest astrophysical objects with (Lewin et al m3 . G allowa^Einf>08). These photospheric 

^ ; apparent surfaces. The supra-nuclear density of the NS in- radius sion (PRE) bursts ide important information 

^ . ner core makes it difficult to construct the theoretical equa- about the NS compactness _ the obserV ed Eddington flux and 

C£ tion of state (EoS), but it can be constrained if the NS masses ±e maximum effective temp erature of the NS surface (e.g . 

O and radii are detenmne d using astrophysical methods (see Ebisuzaki 1987; Dame n"et1l] [19901 Ivan Paradiis et all ^990). 

O . | Lattimer&Prakash| | 2007i for a review). For example, it is pos- ^ first Qne giyes a distance . dependent mass . rad ius relation, 

; sible in principle to constrain the NS parameters from the anal- wh0e ±e second Qne can iye a mass _ radius relation> which is 

ysis of the pulse profiles of rapidly rotating NSs - millisecond inde dent of the distance to the NS . There exist two uncer . 

X-ray pulsars (PouJanen&Gierlmski 2003; Le ahy et al. ||2008). f . _ Rrst nf „„ the chpm{ca] mmnnsirinn nf the at . 
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tainties here. First of all, the chemical composition of the at- 

The apparent NS radius can be obtained from the thermal emis- here of x bursti NSs is not well known . Both; the 

£ , sion of nearby isolated NSs (e.g. |Triimper et al. || 2004|) or NSs Eddington i uminosity and the maximum effective temperature 

in low-mass X-ray binaries (LMXBs) in globular clusters dur- depend Qn k (Lewin et al m3) SecQn(k x _ my bursting NSs 

ing quiescence (e.g.jRutledge et al.||2002a||b|; | Heinke et al. || 2006t eme t tra differ from the b i ackbody at the effective tem- 

| Webb & Barret| | 2007|) . This requires an accurate determination mre The ent tra are close to a diluted b i ackbody 

of the distance to the NS as well as a NS atmosphere model at wkh a color temperature Tq larger than the effective temp erature 

low effective temperatures to predict the emergent spectrum. by the color colTection (or hardness) factor /c . T he reason for 

Additional consttaints can be obtained for the X-ray m i s is Compton s cattering in the upper layers of NS atmospheres 

bursting NSs. These are members of LMXBs with a rela- during the burst (Lo ndon^i][l986; Lapid^i3[l986). The 

tively low accretion rate showing quasi-periodic X-ray bursts yalue of /c depends on the chemical compos i t i on and is a strong 

due to thermonuclear burning of hydrogen and/or helium at f unc ti on of the lum inosity when it approaches the Eddington 

ffigbpttom of the accreted envelope (see e.g. | Lewin et al. | limit dp a vi v e t al.lll99Th . And additionally more uncertainties 

| 1993HStrohmayer & Bildsten| | 2006|) . Sometimes these bursts are ^ because ±Q Eddington i um i nosity depe nds on the gravi- 

snong enough and reach the Eddington luminosity L Edd . In tational redshift of the photosphere ( L ewin ^td]fT993h . which 



Send offprint requests to: V. Suleimanov 



varies during the radius expansion phase. 
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Fast cooling of the NS during the burst with corresponding 
variations of the effective temperature allows us to use the whole 
cooling track to determine NS parameters. The burst spectra are 
usually well described by the blackbody function. The normal- 
ization of the blackbody K is related to the NS radius R and the 
distance as 



where cr e » 0.2(1 + X) cm 2 g 1 is the electron scattering 
(Thomson) opacity, and 



l+z = (l-R s /R) 



-1/2 



(4) 



R$ = 2GM/c is the Schwarzschild radius of the NS. The effec- 
tive temperature T e ff can be expressed via I as 



K = 



<R bh (km) 
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(i) r eff = z 1/4 r, 



Edd, 



(5) 



where /?bt> is the blackbody radius, D\q is the distance in units 
of 10 kpc and z is gravitational redshift at the NS surface. If 
the photospheric radius coincides with the NS radius, the evo- 
lution of K in the cooling tail can only be explained by varia- 
tions of the color-correction. Evolution of f c thus should be re- 
flected in the evolution of K~ 1 ^ 4 with flux dPenninx et alj|1989t 
Ivan Paradiis et al.ll990l :fSuleiman ov et alJ2010l) . Comparing the 
data with theory, it is possible to dete rmine the Eddington flu x 
and the apparent radius R^ — R(l + z) dSuleimanov et alj|2010h . 
For a NS in a globular cluster with known distance these can 
be transformed for the given chemical composition into the con- 
straints on the NS mass and the radius. 

It is clear that this method requires a detailed knowledge 
of how the color-correction varies with luminosity. This behav- 
ior depends mainly on the chemical composition of the NS at- 
mosphere and less on the surface gravity logg. Spectra of X- 
ray bursting NSs have been extensively computed by differ- 
ent groups since more than two decades (|Londonet"al1[ l986: 
Lapi dus et al.lll986t lEbisuzakill 19871; iPaylov et alJll99U Madei 
199UlMadei et alj2004tlMaiczvna et al.ll2005l). However , an ex- 



where the Eddington temperature T^dd is the maximum possible 
effective temperature on the NS surface: 



rpA 8 C 



GMc 



R 2 c 



(1+z). 
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It is defined from the equality of g and the radiation pressure 
acceleration g m d 



4/r 



gmd 
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at the NS surface. Here k v is the opacity per unit mass due to 
free-free and bound-free transitions (i.e. true opacity, which is 
much smaller than electron scattering at L « LnddX an d 47r// v is 
the radiation flux. 

The model atmosphere structure for an X-ray bursting NS 
with given input parameters is described by a set of differential 
equations. The first is the hydrostatic equilibrium equation 



dm 



= g ~ gmd, 



(8) 



dm — —pdz, 



tens ive grid of models is m issing. Recently. iMadei et al.l ([2004) 
and iMaiczvna et al.l (120051) computed a set of models, but for a 
fixed effective temperature and varying logg, while physically 
more motivated would be the models with constant logg and 
varying relative luminosity L/^Edd and, therefore, effective tem- 
perature. 

In the present work we extend previous calculations of NS 
atmosphere models. We consider three realistic NS surface grav- 
ities, various chemical compositions, and a broad range of rela- 
tive luminosities. We compute the spectra as well as the color |Grebenev & Sunvaev lloO 
correction factors as a function of L/Lndd- We show that the de- 
pendence of the blackbody normalization K ~'^ 4 on the observed 
flux can provide information about the ratio of the NS radius to 
the distance, the Eddington flux and the chemical composition 
of the atmosphere. These can be used to constrain the NS pa- 
rameters. 



where P s is a gas pressure, and the column density m is deter- 
mined as 



(9) 



with p denoting the gas density and z the vertical distance. 

The second is the radiation transfer equation with Compton 
scattering taken into c o nsider ati on using the Kom p aneets 



operator (Kompaneetsl ^1957t IZavlin & Shibanovl Il991h 



d 2 fyJy 

drl 



k v + cr. 



(/» - By) ■ 



kT 



k v + CT e TOpC 2 



d I dJy r eff 

X X— \X— - 3Jy + —XJ V 



1 + 



CJy 
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2. Method of calculations 

We compute hot X-ray bursting NS model atmospheres in hydro- 
static and radiative equilibria in plane-parallel approximation. 
Input parameters of the models are the following: chemical com- 
position (especially important is the hydrogen mass fraction X), 
surface gravity 



where x = hv/kT^ is the dimensionless frequency, / v (t v ) * 
1/3 the variable Eddington factor, J v is the mean intensity of 
radiation, B v is the blackbody (Planck) intensity, T is the local 
electron temperature, and C = c 2 h 2 / 2(kT e ff) 3 . The optical depth 
t v is defined as 

dr v — (k v + o~ e ) dm. (11) 
These equations are completed by the energy balance equation 



GM 



-(1+z). 



(2) 



ky (Jy - By) dV ~ <T e 



kT 



and the relative NS luminosity I = L/Lndd- Here R and M are the 
NS radius and the mass, and Lndd is the Eddington luminosity as 
measured at the NS surface: 



4 j JydV-~Y J Xjy 



1 + 



CJy 



dv\ = 0, 



4nGMc 

^Edd = (1 +Z), 

cr. 



(3) 



the ideal gas law 

P* = N m kT, 



(12) 



(13) 
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where iV tot is the number density of all particles, and also by 
the particle and charge conservation equations. Here we as- 
sume local thermodynamic equilibrium (LTE) in our calcula- 
tions, so the number densities of all ionization and excitation 
states of all elements were calculated using Boltzmann and 
Saha equations, but accounting for the pressure ionization effects 
on at omic populations using the occupation prob ability formal- 
ism ([Hummer & Mihalasj 1 1988b as described by Huben v et alj 
(1994j). We take into account electron scattering and free-free 
opacity as well as bound-free transitions for all ions of the 15 
most abundant chemical ele ments (H, He, C, N, O, Ne, Na, Mg, 
Al, Si, S, Ar, Ca, Fe, Ni) (seellbragimov et alj|2003l) using opac- 
ities from Verner & Yakovlevl(ll995l) . 

For solving the abo ve equations, we use d a version of the 
computer code ATLAS dKuruczlll970l ll993), modified to deal 
with high temperatures. The c ode was also modified to ac- 
count for Compton scatteri ng dSuleimanov & Poutanenl f2006: 
ISuleimanov & Wernerll2007l) . 

The course of calculations is as follows. First of all, the 
input parameters of the model atmosphere are defined and a 
starting model using a grey temperature distribution is cal- 
culated. The calculations are performed with a set of 98 
depth points m\ distributed logarithmically in equidistant steps 
from m * 10~ 7 g cirT 2 to m max = 10 6 g crrT 2 . The ap- 
propriate value of m max is such that satisfies the condition 

V T v,b-f,f-f('«max)Tv('«max) > 1 at all frequencies. Here r V]b -f,f-f 
is the optical depth computed with the true opacity only (bound- 
free and free-free transitions, without scattering). Satisfying this 
equation is necessary for the inner boundary condition of the ra- 
diation transfer problem. 

We used the condition of the absence of irradiation flux at 
the outer boundary 



3Jy 

dr v 



(14) 



where h v is the surface variable Eddington factor. The inner 
boundary condition, 



8Jy 

8t v 



8By 

8t v 



is obtained from the diffusion approximation J, 

1/3 XdBy/dTy. 

The boundary conditions along the frequency axis are 



(15) 

B v and H v « 



Jy = By 

at the lower frequency boundary (v n 
kT eff ), and 



(16) 



= 10 14 Hz, hv min « 



X ~ JJy + XJ y 

ox T 







(17) 



10 19 Hz, hv„ 



at the upper frequency boundary (v max 
kT e ff). Condition ( fT6b means that at the lowest energies the 
true opacity dominates over scattering k Y » cr e , and therefore 
J v ~ By. Condition ( fTTI ) means that there is no photon flux along 
the frequency axis at the highest energy. 

For the starting model, all number densities and opacities 
at all depth points and all frequencies are calculated. We use 
300 logarithmically equidistant frequency points in our com- 
putations. The radiation transfer equation ( TTOb is n on-linear 
and i s solved i teratively by the Feautrier method {M ihalas 
19781 see also IZavlin & Shibanovl fl99ll iPavlov et alT Tl991: 
Grebenev & Sunvaevl 120021) . We use the last term of Eq. (TTOb 



in the form xJ' y (l + CJ' v ~ 1 /x 3 ), where J' v ~ l is the mean inten- 
sity from the previous iteration. During the first iteration we take 
J 1 ' 1 = 0. Between iterations we calculate the variable Eddington 
factors f v and h v , using the formal solution of the radiation trans- 
fer equation in three angles at each frequency. We use six iter- 
ations, because usually in the considered models 4-5 iterations 
are sufficient for achieving convergence. 

The solution of the radiative transfer equation dTOb is 
checked for the energy balance equation (TT2l) . together with the 
surface flux condition 



4;r // v (m = 0)a!v = cr S B7; 4 ff = 4nH . 
Jo 

The relative flux error 
H 



s H (m) = 1 - 



JT™ Hy(m)dv 
and the energy balance error as functions of depth 

kT 



(18) 



(19) 



s A (m) = 



ky ( Jy - By) ~ (T g 



TO e C 



4 I JydV — I Xjy 

Jo T Jo 



1 + 



CJy 



dv 



(20) 



are calculated, where H v (m) is the radiation flux at any given 
depth m. This quantity is found from the first moment of the 
radiation transfer equation: 



dfyJy 

<9t v 



Hy 



(21) 



Temperature corrections are then evaluated using three different 
procedures. The first is the integral A-iteration method, modified 
for Compton scattering, based on the energy balance equation 
(fT2l . In this method the temperature correction for a particular 
depth is found from 



AT, 



-e A (m) 



£° [(A v diag - 1)/(1 - a v Ay diag )]ky(dBy/dT)dv 



(22) 



Here a v = cr e /(ky + cr e ), and A v di ag is the diagonal matrix ele- 
ment of the A operator. This procedure is used in the upper at- 
mospheric layers. The second procedure is the Avrett-Rrook flux 
correction, which uses the relative flux error enim) and is per- 
formed in the deep layers. And the third one is the surf ace cor- 
rection , which is based on the emergent flux error. See iKuruczl 
(Il970h for a detailed description of the methods. 

The iteration procedure is repeated until the relative flux er- 
ror is smaller than 1 %, and the relative flux derivative error is 
smaller than 0.01%. As a result of these calculations, we obtain 
a self-consistent X-ray bursting NS model atmosphere, together 
with the emergent spectrum of radiation. 

This code was tested against the code by M adei et al.l d2004h 
on problems of the atmo spheres of X-ray bursting NSs and 
hot DA white dwarfs b y [Suleimanov & Poutanen ( 2006) and 
Suleiman ov et al.l {2006), respectively. See detailed discussion 
in Sect. |5l 



3. Results of atmosphere modeling 

Using the code described above we have calculated an extended 
set of the hot NS atmospheres. The model atmospheres were 
computed for three values of surface gravity logg = 14.0, 14.3 
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6 8 10 12 14 16 
Neutron star radius R (km) 

Fig. 1. The mass-radius relations for several equations of state 
of neutron and strange stars matter are shown by dotted curves. 
The M-R relations for constant \ogg = 14.0,14.3,14.6 are 
shown by solid curves and the lines of constant redshift (z = 
0.18,0.27,0.42 corresponding to the chosen logg for M = 
1.4M ) by dashed curves. The upper-left r egion is excluded by 
constraints from the caus ality requirements dHaensel et al.l2007t 
lLattimer & Pra kash 2007). The hatched horizontal belt marks 
the sprea d of pulsars masses accurately measured in double NS 
binaries (Haensel et al. 2 0071) . 



and 14.6, which cover most of the physically realistic NS equa- 
tions of state for a large range of NS masses (see Fig.Q]). We con- 
sider six chemical compositions: pure hydrogen, pure helium, 
and a solar mixture of the hydrogen and helium with various 
abundances of heavy elements scaled to solar abundances, Z = 
1, 0.3, 0.1 and 0.01 times the solar (Z ), or using the standard 
stellar atmosphere definition [Fe/H] = 0, -0.5, -1 and -2. W e 
used new values of the solar abundances (lAsplund et alj r2009). 
In particular, according to this work the number ratio of he- 
lium to hydrogen is less than was adopted before (n(He)/n(H) 
= 0.0851 instead of 0.0977, a new hydrogen mass fraction X = 
0.7374). For every values of abundance and log g twenty models 
with relative luminosity / = 0.001, 0.003, 0.01, 0.03, 0.05, 0.07, 
0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 
0.98 were calculated. The total number of the computed mod- 
els in the set is 360. Basic properties of the obtained models are 
shown in Figs. [2^41 Most of them are well known from the pre- 
vious investigations (see SectionQ} and we present these figures 
mainly for illustrative purposes. 

The effect of chemical composition on model emergent spec- 
tra and temperature structures is illustrated in Fig. [2] Normalized 
spectra and temperature distributions for models with fixed log g 
= 14.0 and relatively high (/ = 0.95) and low (I = 0.1) luminosi- 
ties for three different chemical compositions (pure hydrogen, 
pure helium and solar abundances) are shown. The dependences 
of the normalized spectra on the photon energy scaled to the ef- 
fective temperatures are presented. 

Atmosphere temperature structures in optically thin layers 
(at electron scattering optical depth r e < 1) are determined by 
the balance between heating of electrons due to down-scattering 
of high energy photons from deep layers and cooling due to 
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Fig. 2. Emergent (unredshifted) spectra and temperature struc- 
tures of model atmospheres with two relative luminosities (I = 
0.95 and 0.1) at fixed surface gravity (logg = 14.0) for various 
chemical compositions: pure hydrogen (dotted curves), pure he- 
lium (dashed curves) and solar mix of hydrogen and helium, and 
subsolar metal abundance Z = 0.3Z o (solid curves). Top panel: 
Spectra normalized to maximum flux versus photon energy in 
units of the effective temperature. Note, that T e ff oc Z'/ 4 (l + 
X)' 1 ^. For clarity, the spectra for I = 0.95 are shifted along 
the ordinate axis by a factor of 10. Bottom panel: Corresponding 
temperature structures (in units of the effective temperature) ver- 
sus Thomson optical depth. 



free-free and bound-free emission. In the upper low-density lay- 
ers, cooling is inefficient and therefore the temperature rises and 
forms a chromosphere-like structure in less luminous models. In 
the high luminosity atmospheres, radiation pressure is signifi- 
cant and the plasma density is low. Due to these conditions, the 
extended high temperature layers stretch up to optically thick 
layers because of the influence of Compton heating. 

Emergent spectra of high temperatu re NS atmospheres 
are close to dilute d blackbody spectra dLondon et all 119861 
lLapidus et al.1l 19861) 

= -^B E (T C = / c r eff ) (23) 

with color correction (or hardness) factor f c (see Section |4] for 
details). To first approximation f c is equal to the ratio of the up- 
per layers (surface) temperature to T e ff . Cooling is more effective 
in pure helium atmospheres, therefore they have lower surface 
temperatures and relatively soft spectra (smaller / c ). 
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Column density (g cm" 2 ) 



Fig. 3. Emergent (unredshifted) spectra (top panel) and temper- 
ature structures (bottom panel) of the model atmospheres with 
four relative luminosities (I = 0.5, 0.1, 0.01 and 0.001) and 
fixed surface gravity (logg = 14.0) for solar hydrogen-helium 
mixture and various abundances of heavy elements: Z = Z 
(solid curves), Z = 0.3Z o (dotted curves), Z = 0.1Z o (dashed 
curves), Z = 0.01Z o (dot-dashed curves). In the top panel, the 
blackbody spectra with effective temperatures are also shown by 
short-dashed curves. 




0.5 1 1 ' 1 ' 1 ' 1 ' 1 

10" 7 10" 5 10" 3 10" 1 10 1 

Column density (g cm" 2 ) 

Fig. 4. Emergent (unredshifted) spectra (top panel) and tem- 
perature structures (bottom panel) for three relative luminosi- 
ties (I = 0.95, 0.5 and 0.1) and fixed chemical abundance (so- 
lar hydrogen-helium mixture and Z = 0.3Z o ) for three differ- 
ent surface gravities: logg = 14.0 (solid curves), 14.3 (dashed 
curves) and 14.6 (dotted curves). For clarity, the spectra for / 
=0.95 and 0.5 are shifted along the ordinate axis by factors 10 +1 ° 
and 10 +0 5 , respectively. 



In the low luminosity atmospheres with heavy elements, iron 
is not completely ionized. Therefore, an absorption edge at 9 
keV due to the bound-free transition from the ground level of 
the H-like Fe xxvi ion arises. This edge reduces the number of 
hard photons, therefore, the upper layers stay at lower tempera- 
tures. The strength of this edge and the surface temperature de- 
pend on the metal abundance as demonstrated in Fig. [3] At the 
lowest luminosity, absorption edges in the 1-3 keV range due 
to other chemical elements arise. It is interesting that emergent 
spectra and temperature structures of the more luminous models 
(/ > 0.5) depend on the heavy element abundances very little, 
because at these high temperatures iron is completely ionized. 

Emergent spectra and temperature structures for three rel- 
ative luminosities and for three different surface gravities are 
shown in Fig. [4] for fixed chemical abundances (solar H/He mix 
and Z = O.3Z ). Again, the luminous models for a given I de- 
pend on log g very little, but for the low luminosity models this 
dependence is significant. The model with lowest gravity has the 
lowest r e jf and the lowest iron ionization. It leads to the largest 
absorption edge and the lowest surface temperature. 



4. Color correction factor 

For a given model atmosphere defined by /, log g and chemical 
composition, we compute also the emergent spectra. These spec- 
tra are then fitted by a diluted blackbody spectrum ( 1231 in some 
energy band, for example, 3-20 keV, corresponding to the very 
often used PCA detector of the RXTE observatory. As in the case 
of the observations, we use the fits with two free parameters, the 
normalization w and the color-correction: 

F E * wB E (f c T etf ). (24) 

This two-parameter approximation gives rather good fits, es- 
pecially for low-temperature atmospheres with heavy elements 
(see Fig. [7]i. The difference between w and f~ 4 is small for lumi- 
nous models. 

Results of the fitting de pend on the use d fitting procedure, 
as was first pointed out by lEbisuzakil d 19871) . Here we use five 
different procedures. In the first one, we minimize the sum 

N 

J](F En - Wl B E Xf c ,iT eff )) 2 , (25) 
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Fig. 5. Dependence of the color correction factors f c on the 
relative luminosity for various NS atmosphere models. The 
f c are obtained using the first fitting procedure. Top panel: 
Dependences for hydrogen and solar H/He mixture with Z = 
0.3Z o models with different surface gravities \ogg = 14.0 (solid 
curves), 14.3 (dotted curves) and 14.6 (dashed curves). For clar- 
ity, the dependences for hydrogen models are shifted up by +0.2. 
Bottom panel: Variation of f c on chemical compositions: pure 
hydrogen (upper curve), pure helium (lowest curve), and solar 
H/He mixture with Z = Z (solid curves), Z = 0.3Z G (dotted 
curves), Z = O.1Z (dashed curves), Z = 0.01Z o (dot-dashed 
curves) for low gravity \ogg = 14.0 models. 



where N is the number of photon energy points in the considered 
energy band. As the energy points in the computed spectra are 
equidistant in logarithm, this procedure is formally equivalent to 
minimizing the integral 



X 



{F E -w x B E (J CtX T eB )f-^. 

£„„„ E 



(26) 



When fitting the data, one fits the photon count flux, not the en- 
ergy flux, therefore, in the second procedure we minimize the 
following sum 



z 

n=\ 



(F En - W2 B En (f ca T eff )) 2 



El 



which is equivalent to minimizing the integral 



X 



(F E -w 2 B E (f ca T eS )) 2 ^. 
E™, E 3 



(27) 
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Fig. 6. Same as bottom panel of Fig. [5] but in a log-scale. 

In the third procedure we suggest to minimize the integral 
JT ™(F E -w 3 B E (f c jT eS )) 2 dE, 
which corresponds to minimizing the sum 

N 

^(F £ „ - w 3 B E Xf, 3 T eff )) 2 E n . 

n=\ 



(29) 



(30) 



We also use a fourth fit with only one free parameter, relating 
w - fTJ, and using the same minimization as for the first proce- 
dure. And finally in the fifth procedure we compute the color cor- 
rection / Cj 5 by dividing the energy where the peak of the model 
flux F E is reached by the peak energy of the blackbod y spectrum 
B E (T e g) as was done by Ma dei et al.l d2004l) and lMaiczvna et all 
(I2005h . 

The obtained color correction factors f c ,i—fcA depend on the 
chosen energy band. We perform our four fitting procedures and 
calculate the corresponding color correctios and dilution factors 
in the 3-20 keV energy band corresponding to the RXTE/PCA 
detector. Due to gravitational redshift a spectrum in the observed 
energy band is radiated in the energy band with blueshifted 
boundaries (3 - 20) x (1 + z) keV. Each NS has its own (a priori 
unknown) gravitational redshift. We calculated redshifts using 
logg and adopting a NS mass equal to 1.4 M (see Eqs. (f2]i and 
©). We obtained R = 14.80, 10.88, 8.16 km and z = 0.18, 0.27, 
0.42 for \ogg = 14.0, 14.3, and 14.6, respectively. Varying the 
mass in the interval M -2M introduces 3, 5, 10% uncertainties 
in 1 +z (see Fig.Q}, which, however, have less than a 0.1% effect 
on the color corrections. The results of the fitting procedures are 
presented in Table [T]and in the online material. 

The obtained results are illustrated in Figs.l5l-fT0l The color- 
correction factors for models with various surface gravities and 
chemical compositions versus relative luminosity are shown in 
Fig. |5] A major feature for all dependences f-L/L^ is a local 
minimum of f c at some intermediate (/ ~ 0.1-0.5) luminosity. 
The color-correction factors are largest at luminosities closest 
to the Eddington luminosities and decrease as the luminosity 
decreases. The reason is well known, it is due to a decreasing 
role of the Compton scattering in comparison with true opac- 
ity (see Section [3). At low luminosity, the color-correction fac- 
tor increases again (see Fig. |6j. At these conditions, Compton 
scattering is not so significant (Suleima nov & Wernerl [2007: 
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Table 1. Color-correction and dilution factors from the blackbody fits to the spectra of hydrogen atmosphere models at log g = 14.0. 
The results for other chemical compositions and gravities are presented in the Online Material. 



Set 1 


X = 1 


Z = lo 


gg= 14.0 T EM 


= 1.644 keV R 


= 14.80 km z = 


0.18 


/ 


T eB (keV) 


fci 


fo,2 


/c,3 


/c.4 


/c.5 




Wo f 4 n 




0.001 


0.292 


1.933 


1.932 


1.934 


1.765 


1.538 


0.742 


0.743 


0.741 


0.003 


0.385 


1.775 


1.774 


1.775 


1.638 


1.448 


0.829 


0.830 


0.829 


0.010 


0.520 


1.648 


1.649 


1.647 


1.543 


1.382 


0.900 


0.900 


0.900 


0.030 


0.684 


1.563 


1.565 


1.560 


1.520 


1.396 


0.942 


0.941 


0.942 


0.050 


0.777 


1.530 


1.533 


1.528 


1.523 


1.429 


0.954 


0.954 


0.955 


0.070 


0.845 


1.511 


1.514 


1.509 


1.517 


1.448 


0.961 


0.961 


0.961 


0.100 


0.924 


1.494 


1.497 


1.491 


1.506 


1.461 


0.966 


0.967 


0.966 


0.150 


1.023 


1.478 


1.480 


1.475 


1.494 


1.469 


0.970 


0.971 


0.970 


0.200 


1.099 


1.470 


1.472 


1.467 


1.487 


1.471 


0.972 


0.973 


0.972 


0.300 


1.217 


1.467 


1.467 


1.465 


1.484 


1.479 


0.975 


0.975 


0.975 


0.400 


1.307 


1.475 


1.474 


1.473 


1.491 


1.491 


0.977 


0.976 


0.976 


0.500 


1.382 


1.491 


1.489 


1.489 


1.506 


1.511 


0.979 


0.978 


0.979 


0.600 


1.447 


1.513 


1.511 


1.511 


1.528 


1.535 


0.980 


0.978 


0.979 


0.700 


1.504 


1.549 


1.547 


1.547 


1.560 


1.571 


0.986 


0.984 


0.985 


0.750 


1.530 


1.567 


1.566 


1.564 


1.578 


1.589 


0.984 


0.984 


0.983 


0.800 


1.555 


1.594 


1.596 


1.591 


1.601 


1.617 


0.992 


0.993 


0.991 


0.850 


1.578 


1.626 


1.630 


1.622 


1.631 


1.647 


0.994 


0.997 


0.993 


0.900 


1.601 


1.668 


1.678 


1.662 


1.672 


1.688 


0.994 


1.000 


0.993 


0.950 


1.623 


1.740 


1.759 


1.731 


1.737 


1.758 


1.004 


1.015 


1.001 


0.980 


1.635 


1.802 


1.832 


1.788 


1.809 


1.815 


0.992 


1.009 


0.989 



Rau ch et al.l l2008>. and the increasing spectral color correction 
is related to the properties of t he free-free opacity as first noted 
by lLondon et all d 19841 119861) . If the relative luminosity / de- 
creases below 0.2, the maximum of the Planck function shifts 
to lower energies. At low photon energies, the free-free opacity 
dominates over electron scattering. In pure hydrogen or helium 
atmospheres, radiation escapes at higher energies easier due to 
the strong, oc E^, energy dependence of the free-free opacity. 
This effect grows as r e (j decreases, and therefore, the spectrum 
becomes relatively harder. In models including heavy elements, 
numerous absorption edges distort this picture, but the common 
opacity properties are the same, the opacity increases as the pho- 
ton energy decreases. 

The computed dependences for model atmospheres with 
heavy elements show relatively deep minima in f c at / ~ 0.1 
(effective temperatures about 1 keV). These minima arise due to 
iron ions, because at those temperatures iron is not completely 
ionized and a significant absorption edge of Fe xxv and xxvi ap- 
pears at ~ 9 keV (Fig. [3}. Due to this discontinuity in the spectra, 
the color correction factors decrease. The depth of the minimum 
depends on the metal abundance. In the case of solar abundance 
(Z = Z©) the minimum is deep and it disappears in the case of Z 
= 0.01 Zq. At the same time the color correction factor at high 
luminosities depends very little on Z. At these luminosities, / c 
depends mainly on a relative hydrogen fraction X, being smaller 
for smaller X, i.e. pure helium models have the smallest color 
correction factors. The dependence of f c on surface gravity is 
also insignificant at high luminosities (see top panel of Fig.|5J- 

Examples of theoretical spectra fitted by the first and the 
fourth (one-parameter) procedures are shown in Fig. [7] For high 
luminosity models, both fits are identical, but at lower luminosi- 
ties the difference is obvious. The fourth procedure gives a bet- 
ter representation of the overall spectral energy distributions, but 
the two-parameters fits (method 1) represent more accurately the 
theoretical spectra in the considered energy band. The difference 
between these two fits is more significant for the model with 
heavy elements due to numerous absorption edges. 



Residuals between theoretical spectra and the two fits, ob- 
tained using the first and the second procedures versus photon 
energy for models with three relative luminosities (/ = 0.1, 0.5 
and 0.95) and two different chemical compositions (X-l and X- 
0.7374, solar composition) are shown in Fig. [8] At photon ener- 
gies >3— 4 keV the fits have an accuracy better than 2-3 % for 
luminous models with heavy elements and for all pure hydro- 
gen models. The residuals for these models increase at energies 
corresponding to the Wien tail (E > 10 keV) due to exponential 
flux decrease. Here the fits overestimate the computed spectra. 
On the contrary, the blackbody fits underestimate the theoretical 
spectra at energies <2-3 keV (see also Fig. [7) and the residuals 
grow here, too. The fits to the spectrum of the low luminosity 
model with heavy elements are worse due to the strong absorp- 
tion edge at m 9 keV. Fitting accuracy of pure helium model 
atmosphere spectra is similar to pure hydrogen model spectra. 
We can conclude that at high luminosities (/ > 0.5) the diluted 
blackbody spectra are rather good representations of the theoret- 
ical NS atmosphere model spectra within the RXTE/PCA energy 
band, with deviations not exceeding 2%. 

Absolute and relative comparisons of color corrections / c 
and corresponding dilution factors w c which were obtained using 
all four procedures are shown in Figs. [9] and [10] The 2-parameter 
fitting procedures 1-3 give almost identical results differing in 
/ c by not more than 2%. The 1-parameter procedures 4-5 give 
also very similar / c in the luminosity range / > 0.2. The differ- 
ences between 1- and 2-parameter procedures become signifi- 
cant at lower luminosities, because these spectra strongly devi- 
ate from a diluted blackbody (these deviations are significantly 
larger for models w ith metals). The fifth proc edure employed by 
Mad ei et al.1 (12004]) and lMaiczvna et afl (120051) gives values of / c 
at I < 0.2 significantly lower than procedures 1-3. 

5. Comparison with previous calculations 

There is an extensive literature on NS atmosph e res. T he first 
detailed models are presented by lLondon et al.l (119841) . They 
considered hydrogen-helium-iron atmospheres with the relative 
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1 10 

E (keV) 

Fig. 7. Examples NS atmosphere spectra for high (7=0.95) 
and low (7=0.001) luminosity and low gravity (logg = 14.0). 
Theoretical model spectra are shown by the solid curves. The di- 
luted blackbody fits with the first fit procedure are shown by the 
dashed curves and the one-parameter fits (fourth fit procedure) 
are shown by the dotted curves. The vertical dotted lines show 
boundaries of the energy band where the fitting procedure was 
performed. Top panel is for hydrogen atmosphere and the bottom 
panel is for the solar H/He mixture with Z = 0.3Z o . 



number density ratio He/H=0.1/1 and various iron abundances 
(Fe/H=3.4xl0~ 5 =solar, 0.1 solar and 10~ 5 of solar) and ac- 
counted for Compton scattering using the Fokker-Planck ap- 
proximation (Kompaneets equation). They showed that/ c (com- 
puted using 1 -parameter procedure 4) is large at both high and 
low ends of relative luminosities. If the metal abundance is low, 
f c shows a broad minimum at / w 0.2 of about 1.4, while 
the models with solar Fe abundance also show a dip in f c at 
/ 0.1. The presen ce of the dip is also discussed in detail by 
Lapidusetal . (1986). These results are consistent with our cal- 
culations at qualitative level. 

At luminosities close to Eddington. lPavlov et al. | (Eillb de- 
rived a simple approximate formula for the ratio of surface tem- 
perature (which is not far from the color temperature) to T e s 

For our luminous models presented on Fig. [2] (7 = 0.95, 
logg = 14.0) this equation gives r surf /r eff = 1.518 (X=0), 1.569 
(X=0.7374), and 1.582 (X=l). In our models we have 1.534, 
1.609 and 1.641, respectively. The relative deviations are 1.1, 
2.5 and 3.7%. 




10 

E (keV) 



Fig. 8. Relative errors of the computed spectra fitted by the first 
(solid curves) and the second (dashed curves) fit procedures vs. 
photon energy for hydrogen {top panel) and solar H/He mixture 
with Z = 0.3Z o {bottom panel) low gravity (logg =14.0) mod- 
els. Corresponding relative luminosities and effective tempera- 
tures are given at the curves. The vertical dotted line shows the 
low boundary of the energy band where fitting procedures were 
performed. For clarity, models with / = 0.5 and 0.1 were shifted 
by +0. 1 and +0.2, respectively. 

We also compared our results for low temperature (T e ff=l-3 
MK, when Compton scattering is insignificant), pur e H and pure 
He atm osphere models to previous calcu lations b y Zavlin et aL 
(1 19961) and found a good agreement (ISuleimanov & W erner 
2007). Comparison of our LTE models with heavy elements 
to the non-LTE m odels computed usin g the TMAP cod e 
dWerner et al.l l2003) has been performed bv lRauch et~ al. (2008). 
It confirmed the reduction of the iron absorption edge at 9 keV 
for models with T e ff 1 keV due to iron over-ionization, as was 
found before bv lLondon et al.l (Il986l) . However, the non-LTE ef- 
fects reduce the edge strength by one third only, much less than 
was claimed by London et al.l (1 1 9861) . who considered only one 
level for Fe xxvi. 

Recently, iMaiczvna et alj d2005l) have computed a set of 
atmosphere models for a small set of effective temperatures 
and various \ogg = 12.9,..., 15.0. They considered hydrogen- 
helium-iron atmospheres with the relative number density ra- 
tios Fe/He/H=3.7xlO~ 5 /0-ll/l corresponding to the 'old' so- 
lar abundances with X=0.693. They assumed LTE, ignored 
pressure ionization, and employed the integral formulation of 
Compton scattering using isotropic redistribution function. The 
radiative transfer equation was solved simultaneously with ra- 
diative equilibrium equation by the iterative complete lineariza- 
tion method. For log = 14.6 (r Edd = 2.42 keV for M = 
1.4M Q ) they obtained f c = 1.33,1.26,1.40,1.54 for T e $ = 
(1, 1.5, 2, 2.5) x 10 7 K which correspond to relative luminosities 
/ = 0.016, 0.081, 0.26, 0.63. Our calculations (applying the fifth 
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Fig. 9. Top panel: Dependence of the color correction factors 
on the relative luminosity for the low surface gravity (logg = 
14.0) calculated by the five fitting procedures. The results ob- 
tained with the first procedure are shown by the solid curve, the 
second procedure - by the dashed curve, the third - by the dot- 
dashed curve, the fourth - by the dotted curve, and the fifth - 
by circles. The lower curves are for the solar mixture of hydro- 
gen and helium and Z = 0.3Z o . The upper curves correspond 
to pure hydrogen models and are shifted up by +0.2. Bottom 
panel: Ratio of the color correction factors obtained using the 
second (dashed curve), third (solid curve), fourth (dotted), and 
fifth (circles) procedures to the color correction factor from the 
first procedure. The curves corresponding to hydrogen models 
are shifted up by +0.2. 
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Fig. 10. Same as in Fig. [9] but for the dilution factor. The results 
are shown for the procedures The curves corresponding to 
hydrogen models are shifted up for clarity. 



It is also necessary to notice that obtaining f c for low relative 
luminosity and low gra vity (e.g. logj? = 14.0) from the results 
presented in Table 1 of Maiczvna et al.l (120051) (that have mini- 
mum I = 0.06) is impossible. One cannot scale the results ob- 
tained for logg = 15.0, because the opacity due to free-free and 
bound-free transitions depends on temperature and density in a 
different way, resulting for the same I in smaller ionization and 
stronger absorption edges at higher logg (see Fig. @J, affecting 
thus the color temperature estimation. 



fit procedure used by Maicz vna et alj |2005) for the same logg, 
I, and for X = 0.74, Z = Zq (see Set 6 in Table 1 of the Online 
Material) give / c , 5 = 1.19, 1.27, 1.39, 1.52. We see that at high 
luminosities, the difference is about 1%, while at the lowest tem- 
perature our color correction is significantly smaller. This can be 
explained by the difference in the chemical abundance, since we 
have included among other elements also Ne, Mg, Al, and S, 
which produce edges in the 1-3 keV energy range strongly af- 
fecting f c . Note that our pure hydrogen atmosphere model which 
lacks the edges gives / c .g = 1.39 for / = 0.016. In general, we 
find the agreement rather satisfactory. 

We remark, however, that the color-corrections obtained 
from the fifth procedure should not be compared to the data, be- 
cause the color temperatures of the time-resolved spectra from 
X-ray bursts are computed by fitting the actual data in a specific 
energy interval (e.g. 3-20 keV) with the diluted blackbody func- 
tion with arbitrary normalization. Thus, it is more appropriate to 
use color-correction factors obtained by procedures 1-3 (which 
give almost identical results). 



6. Application to observational data 

6.1. Method 

We discuss now how our theoretical calculations can be used for 
comparison with the observational data. Observed spectra of X- 
ray bursting NSs are usually fitted by the blackbody (see e.g. 
Gallo way etai]l2008h with two free parameters: a color temper- 
ature 7\,b and a normalization K = (R\,b [km]/Dio) 2 (here 
is the apparent blackbody radius measured in km and Dio is a 
distance to the source in units of 10 kpc). A total observed flux 
F - R^ b os,BT* b / D 2 can be also expressed using basic NS pa- 
rameters F = # 2 o- SB 7^ ff (l + zT 2 ID 2 dLewinet al.lll993l) . Using 
the obvious relation 7bb = f c T s «(X + z) _I we get 

Rco=R(l+z)= Rbbfl (32) 

where Roa is the apparent NS radius at infinity. We can rewrite 
Eq. (l32l as 

f c =K- 1/4 A, (33) 
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Fig. 11. The dependence K~ l ^ 4 -F as observed during the cool- 
ing track of the long burst from 4U 1724-307 on November 8, 
1996 (stars). The theoretical f c -l dependence is shown by the 
dashed curve (right and upper axes) and the best-fit relation 
(solid curve). 
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Fig. 12. Constraints on M and R from various observables. The 
solid curve gives the relation obtained from the Eddington tem- 
perature d36l l using Equations ( 1371 ) and J38i . Dotted curves are 
for the Eddington flux d39l . dashed curves are for A=const, 
shown for two different distances. For the distance below the up- 
per limit (l4Tb there are solutions at the interception of the curves. 
If the assumed distance is too large, there are no solutions (the 
curves shown by thin lines do not cross). 



where 



^oo [km] 

t Dm 



-1/2 



R[km](l+z) 



D 



10 



•1/2 



(34) 



is a constant. 

lb isuzakil (1 1987b suggested to fit cooling tracks presented 
in the form F/Fzdd-Tbb by the theoretical relation /-/ c r e ff. 
However, the flux-temperature dependence is dominated by the 
approximate F oc T^ b relation and therefore it is more appealing 
to emphasi ze the deviations from the constant apparent radius. 
In addition. lEbisuzakH (11987) fixed a priori the value of F E dd to 
the flux at the touchdown point, which corresponds to minimum 
Rbb and maximum T\,b (and flux). Such a restriction is not sup- 
ported theoretically as the maximum flux observed from PRE 
X-ray bursts can exceed the flux when the photo sphere is at the 
stellar surface by a factor 1 +z dLewin et alj!993l) . A small appar- 
ent radius at the touchdown does not necessarily mean that the 
photosphere actually coincides with the NS surface, because at 
luminosities very close to Eddington the col or correction can be 
rather large dPavlov et alJll99lh . Recently, ISteiner et all d201dh 
also argued in favor of this interpretation on the basis of the in- 
consistencies between the observables in some PRE bursts. 

Because the evolution of K~ 1 ^ 4 at sub-Eddington luminosi- 
ties (at late burs t stages) reflects th e evolution of the color 
correction factor ( Penninx et alJI 19891 : Ivan Paradjisetal][l990; 
ISuleimanov et alJl2010h . we propose here to fit the observed de- 
pendence K~ l t 4 -F by the theoretical relations / c -(/ = L/L E dd = 
F/FEdd)- The two free parameters of the fit are A and the (ob- 
served) Eddington flux 



F E dd 



^ 2 



GMc 



AnD 2 



D 2 



(i+zy 



(35) 



(see Fig.rTTTi. Because the evolution of f c is strongest near the 
Eddington flux, the PRE X-ray bursts are most suitable for the 
analysis. We would like to emphasize that the Eddington flux 
should be obtained from the fit to evolution of K~ 1 ^ in a broad 
range of luminosities. 



Combining A and F E dd we can obtain the effective 
(Eddington) temperature corresponding to the Eddington flux on 
the NS surface corrected for the gravitational redshift: 



?Edd,c 



-^V /4 — = 6.4 x 10 9 F l J 4 A' 1 K. (36) 
cr SB crJ l+z Edd 



This quantity is independent of the (uncertain) distance to the 
source and can be used to express the NS radius through the 
observables and the compactness u = Rs/R — 1 - (1 + zY 2 : 



R = 



u(l -m) 3/2 , 



2o- e o- SB r Eddi0o 

and the mass is then found via (see solid curve in Fig.[T2b 

M R 
Mq ~ 2.95 km 



(37) 



(38) 



If the distance D is known (for example, if the source is in a 
globular cluster), one can get additional constraints on M and R 
(for a given chemical composition) from A and ^Edd separately. 
From the Eddington flux estimate we have (see dotted curves in 

Fig.o 



2cr e D F E dd _i ._i/2 
R = 5 u (1 - u) 1 

c j 

= 14.138 km (1 + X) D 2 W F Edd ,-7 u x (1 - i/r 1/2 , 



(39) 



where F E dd-i is the Eddington flux in units 1CT 7 erg crrT 2 s _1 
and the mass is found using Equation (l38l . A measurement of A 
gives another constraint: 



(40) 



R = R m Vl — u = Dio A~ 2 Vl — u km. 



Combining with the parametric expression for the mass (1381 1. we 
get the third relation between M and R shown by the dashed 
curves in Fig. [12] 
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All three curves cross in one or two points (see Fig.[T2l) if the 
distance satisfies the following relation: 

0.0177 

Dm < rr-rrr ■ 41 

(1 + X) A 2 F^dd,-! 

In the opposite case, there is no physical solution for M and R 
for given observables. 

6.2. Comparison to the data 

Although the analysis of the X-ray burst data is outside of the 
scope of the present paper, we give now a few examples of the 
observed bursts that can be interpreted using our approach. We 
also discuss a few cases where the data do not follow the theory 
an d discuss the possible reasons for that. 

iPenninx et al.l d!989l) analyzed two long-duration (>100 s) 
PRE bursts observed by ME/EXOSAT in 1984 and 1986 from 
4U 1608-52 in its hard state at a rather low persistent flux of 
(l-2)xl0~ 9 erg crrT 2 s _1 in 2-20 keV band. The evolution of 
7bb^~^ 4 (which is proportional to K~ 1 ^ 4 ) with \ogF shown in 
their Fig. 7 is almost identical to our models presented in Fig. [6] 
Observations by LAC/Ginga of a long PRE burst from the atoll 
source 4U2129+11 during its island (hard) state at a low per- 
sis tent flux level of ~ 0.5 X 10~ 9 erg cirT 2 s _1 are presented 
by Ivan Paradiis et al.l (11990b . The behavior of T^F ' A at fluxes 
above 30% of the peak (touchdown) flux shown in their Fig. 10 is 
very similar to that in our Fig. [5] A long PRE burst was observed 
by RXTE on November 8, 1996 from 4U 1724-307 in the hard 
state at a l ow accretion rate w i th persistent flux of 1 .2 x 10~ 9 erg 
cirr 2 s 1 (Mol kovet al.ll2 000; Galloway et al.ll2008l) . The evolu- 
tion of K~ 1 / 4 with flux during th e cooling tail (shown in Fig.fTTl 
see also Sulei manov et al.l2010l) is almost identical to that of the 
burst from 4U2129+11. For both objects the data at high fluxes 
are well described by the theory. We note that in both cases the 
position of the touchdown point in K~ x ^-F diagram is not con- 
sistent with extrapolation of the data from intermediate fluxes, 
implying that the Eddington flux (at the NS surface) is actually 
smaller by about 15% than the touchdown flux. This is consis- 
tent with the predictions that the maximum flux during a PRE 
burst (corresponding to the expanded photos phere) can be 1 + z 
times larger than the surface Eddington flux ( Lewin et al ] l!993h . 

However, not all bursts show an evolution consistent with 
theory. For example, two short PRE bursts from 4U 1724-307 
detected by R XTE on Feb 23 and May 22, 2004 during the 
high/soft state dGallowav et al.ll2008h s how a nearly constant ap - 
parent area K during the cooling stage (Suleimanov et al. 2010). 
Four short PRE bursts from 4U 1608-52 at a high persistent flux 
level of (3-6)xl0~ 9 erg cirr 2 s _1 also demonstrate nearly con- 
stant area. This behavior is inconsistent with theoretical expec- 
tations. Therefore, we warn against using this kind of bursts for 
determining NS p arameters such as their masses and radii (cf. 
iGiiver et alll2010h . 

The most plausible reason for deviations of the data from 
theory is an influence of accretion on the emergent spectrum. 
The spectra and variability of the accreting NSs in the soft state 
are consistent with the presence of an op tically thick bound- 
ary/spreading layer at the N S surface dGilfanov et al.l 12003; 
iRevnivtsev & Gilfanovl 12006). The observed persistent spectra 
are well reproduced in the framewo r k of the spreading layer 
mode l dlnogamov & Sun vaev 19991; ISuleimanov & Poutanenl 
2006). During the short bursts, accretion most probably persists 
and therefore photons escape from the atmosphere of a rapidly 
rotating spreading layer. The combined influence of the radia- 
tion pressure and centrifugal force cause a strong reduction of 



the effective gravity (Inogamov & Sunvaevll 19991) which causes 
the spectrum of the escapin g radiation to have a high va lue of 
/ c similar to that at I ~ 1 ( ISuleimanov & Poutan en 2006). This 
effect may be responsible for nearly constant apparent areas in 
spite of the changes in the burst luminosity. 

7. Conclusions 

In this work we present an extended set of the X-ray bursting NS 
model atmospheres in the luminosity range 0.001 - 0.98 Z,Edd for 
three values of surface gravities (logg = 14.0, 14.3, and 14.6), 
and six atmosphere chemical compositions (pure H and He, so- 
lar H/He mixture with Z = 1, 0.3, 0.1 and 0.01 Z Q ). Altogether 
360 models were computed. The calculated models are in a good 
agreement with models presented by other authors at the same 
parameters. 

The spectra of all models are redshifted and fitted by diluted 
blackbody spectra in the RXTE/PCA energy band 3-20 keV. We 
use five different fitting procedures. At relatively large luminosi- 
ties (L > 0.3L£dd) all procedures give almost identical (differ- 
ence less than 2 %) color correction f c - T c /T e ff. The color cor- 
rection factor strongly depends on the relative luminosity and 
the chemical composition and less significantly on the surface 
gravity. For luminous models, color corrections depend mainly 
on the hydrogen mass fraction X. Largest f c are exhibited by 
pure hydrogen models (X = 1). The local minimum in the f c — 
L/Lzdd dependence appears at about 10% Eddington luminosity 
for models with metals due to an appreciable absorption edge at 
9 keV arising from bound-free transitions in hydrogen-like iron. 
The depth of this minimum depends on the metal abundances 
and it disappears in models with Z = 0.01 Z Q . 

The theoretical emergent spectra and the color corrections 
with corresponding normalization are available onlineQ These 
models can be used for interpretation of the X-ray bursting NS 
spectra. We describe in detail the method for the evaluation of 
NS masses and radii using the observed K~ x ^—F dependence 
during cooling tails of PRE bursts. We also discuss the data that 
seem to follow the theory and also the bursts that contradict the 
theory. We argue that only the bursts that clearly follow theoret- 
ical models should be considered for further analysis with the 
aim to determine NS masses and radii. 
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Table 1. Color-correction and dilution factors from the blackbody fits to the atmosphere model spectra. 



1 


T eS (keV) 


Ai 


/c,2 


/c,3 


/c,4 


/c,5 


H<1 


W 2fc.2 


W ifc.3 


Set 1 X = 1 


Z = 


logg = 


14.0 r Edd = 1.644 keV 


R = 14.80 km 


z = 0.18 


0.001 


0.292 


1.933 


1.932 


1.934 


1.765 


1.538 


0.742 


0.743 


0.741 


0.003 


0.385 


1.775 


1.774 


1.775 


1.638 


1.448 


0.829 


0.830 


0.829 


0.010 


0.520 


1.648 


1.649 


1.647 


1.543 


1.382 


0.900 


0.900 


0.900 


0.030 


0.684 


1.563 


1.565 


1.560 


1.520 


1.396 


0.942 


0.941 


0.942 


0.050 


0.777 


1.530 


1.533 


1.528 


1.523 


1.429 


0.954 


0.954 


0.955 


0.070 


0.845 


1.511 


1.514 


1.509 


1.517 


1.448 


0.961 


0.961 


0.961 


0.100 


0.924 


1.494 


1.497 


1.491 


1.506 


1.461 


0.966 


0.967 


0.966 


0.150 


1.023 


1.478 


1.480 


1.475 


1.494 


1.469 


0.970 


0.971 


0.970 


0.200 


1.099 


1.470 


1.472 


1.467 


1.487 


1.471 


0.972 


0.973 


0.972 


0.300 


1.217 


1.467 


1.467 


1.465 


1.484 


1.479 


0.975 


0.975 


0.975 


4.00 


1 ^07 


1.475 


1.474 




1 4.Q1 


1 4.Q1 


Q77 
\j.y i i 


Q7fS 


0.976 


0.500 


1.382 


1.491 


1.489 


1.489 


1.506 


1.511 


0.979 


0.978 


0.979 


0.600 


1.447 


1.513 


1.511 


1.511 


1.528 


1.535 


0.980 


0.978 


0.979 


0.700 


1.504 


1.549 


1.547 


1.547 


1.560 


1.571 


0.986 


0.984 


0.985 


0.750 


1.530 


1.567 


1.566 


1.564 


1.578 


1.589 


0.984 


0.984 


0.983 


0.800 


1.555 


1.594 


1.596 


1.591 


1.601 


1.617 


0.992 


0.993 


0.991 


0.850 


1.578 


1.626 


1.630 


1.622 


1.631 


1.647 


0.994 


0.997 


0.993 


0.900 


1.601 


1.668 


1.678 


1.662 


1.672 


1.688 


0.994 


1.000 


0.993 


0.950 


1.623 


1.740 


1.759 


1.731 


1.737 


1.758 


1.004 


1.015 


1.001 


0.980 


1.635 


1.802 


1.832 


1.788 


1.809 


1.815 


0.992 


1.009 


0.989 



Set 2 X = 1 Z = logg = 14.3 T Edd = 1.954 keV R = 10.88 km z = 0.27 



0.001 


0.347 


1.889 


1.888 


1.890 


1.735 


1.531 


0.777 


0.779 


0.776 


0.003 


0.457 


1.737 


1.736 


1.737 


1.611 


1.443 


0.855 


0.857 


0.855 


0.010 


0.618 


1.615 


1.615 


1.614 


1.525 


1.382 


0.916 


0.916 


0.917 


0.030 


0.813 


1.532 


1.534 


1.530 


1.513 


1.403 


0.951 


0.951 


0.952 


0.050 


0.924 


1.501 


1.504 


1.499 


1.505 


1.430 


0.961 


0.961 


0.962 


0.070 


1.005 


1.484 


1.486 


1.481 


1.496 


1.442 


0.966 


0.967 


0.966 


0.100 


1.099 


1.467 


1.470 


1.464 


1.482 


1.449 


0.970 


0.971 


0.970 


0.150 


1.216 


1.453 


1.455 


1.451 


1.469 


1.452 


0.973 


0.974 


0.973 


0.200 


1.306 


1.448 


1.448 


1.445 


1.464 


1.454 


0.975 


0.975 


0.974 


0.300 


1.446 


1.449 


1.447 


1.447 


1.464 


1.463 


0.977 


0.975 


0.976 


0.400 


1.554 


1.460 


1.456 


1.459 


1.475 


1.480 


0.979 


0.976 


0.978 


0.500 


1.643 


1.479 


1.475 


1.478 


1.493 


1.500 


0.981 


0.978 


0.981 


0.600 


1.719 


1.503 


1.499 


1.502 


1.518 


1.526 


0.980 


0.978 


0.980 


0.700 


1.787 


1.539 


1.537 


1.537 


1.548 


1.562 


0.987 


0.986 


0.987 


0.750 


1.818 


1.562 


1.561 


1.560 


1.572 


1.584 


0.987 


0.986 


0.986 


0.800 


1.848 


1.588 


1.590 


1.584 


1.594 


1.610 


0.992 


0.993 


0.991 


0.850 


1.876 


1.621 


1.627 


1.616 


1.624 


1.641 


0.997 


1.000 


0.995 


0.900 


1.903 


1.665 


1.676 


1.658 


1.667 


1.683 


0.998 


1.005 


0.996 


0.950 


1.929 


1.738 


1.761 


1.728 


1.731 


1.752 


1.007 


1.021 


1.004 


0.980 


1.944 


1.807 


1.840 


1.794 


1.788 


1.817 


1.023 


1.043 


1.020 



Set 3 X=l Z = logg = 14.6 T Edd = 2.322 keV R = 8.16 km z = 0.42 



0.001 


0.413 


1.851 


1.850 


1.851 


1.715 


1.525 


0.804 


0.805 


0.804 


0.003 


0.543 


1.703 


1.703 


1.702 


1.595 


1.438 


0.876 


0.876 


0.876 


0.010 


0.734 


1.584 


1.585 


1.582 


1.512 


1.383 


0.930 


0.929 


0.931 


0.030 


0.966 


1.504 


1.507 


1.501 


1.490 


1.407 


0.959 


0.958 


0.959 


0.050 


1.098 


1.474 


1.478 


1.471 


1.479 


1.425 


0.967 


0.967 


0.967 


0.070 


1.194 


1.457 


1.461 


1.454 


1.468 


1.431 


0.970 


0.971 


0.970 


0.100 


1.306 


1.443 


1.446 


1.440 


1.456 


1.434 


0.973 


0.974 


0.973 


0.150 


1.445 


1.431 


1.433 


1.429 


1.446 


1.435 


0.975 


0.976 


0.975 


0.200 


1.553 


1.428 


1.428 


1.426 


1.443 


1.438 


0.976 


0.976 


0.976 


0.300 


1.718 


1.433 


1.431 


1.432 


1.448 


1.449 


0.978 


0.977 


0.978 


0.400 


1.847 


1.449 


1.445 


1.448 


1.462 


1.468 


0.981 


0.978 


0.980 


0.500 


1.952 


1.470 


1.466 


1.469 


1.484 


1.492 


0.981 


0.978 


0.980 


0.600 


2.043 


1.498 


1.495 


1.496 


1.507 


1.520 


0.987 


0.985 


0.986 


0.700 


2.124 


1.535 


1.534 


1.533 


1.544 


1.556 


0.987 


0.986 


0.987 


0.750 


2.161 


1.558 


1.559 


1.555 


1.561 


1.579 


0.996 


0.996 


0.995 


0.800 


2.196 


1.584 


1.589 


1.579 


1.587 


1.604 


0.995 


0.998 


0.994 


0.850 


2.229 


1.619 


1.627 


1.613 


1.619 


1.637 


1.000 


1.005 


0.999 
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Table 1. Continued. 
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y e ff (Kevj 


/c,l 


Jc,2 


Jc,3 


JcA 


Jc,5 


W l/c,l 




•1, f 4 


0.900 


2.261 


1.664 


1.679 


1.656 


1.662 


1.680 


1.002 


1.012 


1.000 


0.950 


2.292 


1.733 


1.760 


1.721 


1.724 


1.744 


1.011 


1.028 


1.008 


0.980 


2.310 


1.810 


1.852 


1.794 


1.798 


1.816 


1.014 


1.039 


1.010 


Set 4 


X = 0.74 


Z = Z 




= 14.0 


?Edd 


1.703 keV 


R = 


14.80 km 


z = 0.18 


0.001 


0.303 


2.157 


2.148 


2.165 


1 . 1 87 


0.976 


0.096 


0.097 


0.096 


0.003 


0.398 


1.740 


1.724 


1.750 


1.130 


1.176 


0.273 


0.277 


0.270 


0.010 


0.538 


1.448 


1.431 


1.459 


1.156 


1.240 


0.638 


0.648 


0.631 


0.030 


0.709 


1.368 


1.365 


1.366 


1.266 


1.203 


0.890 


0.891 


0.891 


0.050 


0.805 


1.347 


1.361 


1.336 


1.326 


1.129 


0.958 


0.956 


0.962 


0.070 


0.876 


1.328 


1.352 


1.310 


1.326 


1.264 


0.989 


0.988 


0.993 


0.100 


0.957 


1.310 


1.340 


1.290 


1.309 


1.257 


1.000 


1.003 


1.004 


0.150 


1.060 


1.317 


1.342 


1.300 


1.322 


1.317 


0.986 


0.992 


0.987 


0.200 


1.139 


1.341 


1.356 


1.331 


1.354 


1.348 


0.975 


0.979 


0.974 


0.300 


1.260 


1.392 


1.395 


1.390 


1.405 


1.412 


0.979 


0.981 


0.979 


0.400 


1.354 


1.428 


1.423 


1.428 


1.446 


1.447 


0.973 


0.971 


0.973 


0.500 


1.432 


1.457 


1.451 


1.457 


1.473 


1.478 


0.976 


0.973 


0.976 


0.600 


1.498 


1.484 


1.480 


1.483 


1.497 


1.508 


0.982 


0.979 


0.981 


U. /UU 


I. J J / 


JO 10 


JO Jo 


1 <n 
l.M j 


JO 23 


lo4z 


u.ys / 




u.yso 


0.750 


1.584 


1.535 


1.534 


1.532 


1.543 


1.561 


0.989 


0.988 


0.988 


0.800 


1.610 


1.555 


1.556 


1.552 


1.562 


1.582 


0.990 


0.991 


0.990 


0.850 


1.635 


1.588 


1.592 


1.583 


1.595 


1.613 


0.990 


0.992 


0.989 


0.900 


1.658 


1.635 


1.643 


1.629 


1.634 


1.658 


1.001 


1.006 


1.000 


0.950 


1.681 


1.700 


1.716 


1.692 


1.699 


1.720 


1.002 


1.011 


1.000 


0.980 


1.694 


1.750 


1.769 


1.741 


1.731 


1.767 


1.024 


1.036 


1.022 



Set 5 X = 0.74 Z = Z log g = 14.3 T Edd = 2.024 keV R = 10.88 km z = 0.27 



0.001 


0.360 


2.001 


1.973 


2.015 


1.265 


1.298 


0.208 


0.214 


0.205 


0.003 


0.474 


1.651 


1.622 


1.667 


1.204 


1.272 


0.445 


0.459 


0.438 


0.010 


0.640 


1.427 


1.409 


1.433 


1.221 


1.073 


0.771 


0.781 


0.767 


0.030 


0.842 


1.329 


1.347 


1.314 


1.305 


1.081 


0.964 


0.959 


0.971 


0.050 


0.957 


1.287 


1.319 


1.265 


1.289 


1.231 


1.004 


1.002 


1.012 


0.070 


1.041 


1.269 


1.304 


1.248 


1.270 


1.269 


1.006 


1.006 


1.010 


0.100 


1.138 


1.277 


1.305 


1.261 


1.280 


1.283 


0.986 


0.990 


0.987 


0.150 


1.259 


1.313 


1.326 


1.306 


1.328 


1.323 


0.966 


0.970 


0.965 


0.200 


1.353 


1.351 


1.353 


1.351 


1.365 


1.368 


0.976 


0.976 


0.976 


0.300 


1.498 


1.400 


1.393 


1.404 


1.413 


1.418 


0.980 


0.976 


0.981 


0.400 


1.609 


1.428 


1.419 


1.429 


1.444 


1.447 


0.976 


0.971 


0.977 


0.500 


1.702 


1.452 


1.445 


1.451 


1.465 


1.477 


0.981 


0.977 


0.980 


0.600 


1.781 


1.477 


1.472 


1.475 


1.487 


1.506 


0.986 


0.983 


0.985 


0.700 


1.851 


1.507 


1.505 


1.505 


1.518 


1.537 


0.985 


0.983 


0.984 


0.750 


1.883 


1.529 


1.528 


1.526 


1.536 


1.558 


0.991 


0.990 


0.990 


0.800 


1.914 


1.556 


1.556 


1.552 


1.561 


1.582 


0.992 


0.992 


0.991 


0.850 


1.943 


1.586 


1.591 


1.582 


1.589 


1.612 


0.996 


0.999 


0.995 


0.900 


1.971 


1.629 


1.638 


1.622 


1.629 


1.649 


1.000 


1.005 


0.998 


0.950 


1.998 


1.692 


1.710 


1.683 


1.691 


1.707 


1.001 


1.013 


0.999 


0.980 


2.013 


1.750 


1.774 


1.740 


1.729 


1.762 


1.026 


1.041 


1.023 



Set 6 X = 0.74 Z = Z B log g = 14.6 T Edd = 2.405 keV R = 8.16 km z = 0.42 



0.001 


0.428 


1.946 


1.932 


1.951 


1.357 


1.101 


0.327 


0.332 


0.325 


0.003 


0.563 


1.636 


1.627 


1.634 


1.281 


1.203 


0.572 


0.577 


0.573 


0.010 


0.761 


1.397 


1.412 


1.384 


1.284 


1.172 


0.871 


0.862 


0.880 


0.030 


1.001 


1.248 


1.289 


1.225 


1.264 


1.187 


1.027 


1.014 


1.041 


0.050 


1.137 


1.223 


1.260 


1.202 


1.227 


1.250 


1.014 


1.010 


1.024 


0.070 


1.237 


1.239 


1.266 


1.226 


1.238 


1.250 


0.985 


0.987 


0.988 


0.100 


1.352 


1.275 


1.288 


1.272 


1.291 


1.289 


0.961 


0.963 


0.961 


0.150 


1.497 


1.328 


1.325 


1.333 


1.354 


1.331 


0.953 


0.952 


0.953 


0.200 


1.608 


1.363 


1.353 


1.369 


1.388 


1.368 


0.959 


0.955 


0.960 


0.300 


1.780 


1.404 


1.395 


1.407 


1.416 


1.418 


0.982 


0.977 


0.982 


0.400 


1.913 


1.424 


1.418 


1.424 


1.436 


1.447 


0.982 


0.978 


0.982 


0.500 


2.022 


1.444 


1.441 


1.443 


1.455 


1.475 


0.985 


0.982 


0.984 


0.600 


2.117 


1.470 


1.468 


1.468 


1.478 


1.500 


0.988 


0.987 


0.988 


0.700 


2.200 


1.504 


1.503 


1.501 


1.511 


1.531 


0.990 


0.990 


0.990 


0.750 


2.238 


1.526 


1.528 


1.522 


1.529 


1.552 


0.996 


0.996 


0.995 
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Table 1. Continued. 



/ 


1 eff t Ke v ) 


/ci 


/c,2 




JcA 


/c,5 




W 2/ c , 2 


f4 


0.800 


2.275 


1.557 


1.560 


1.553 


1.553 


1.579 


1.005 


1.007 


1.004 


0.850 


2.309 


1.584 


1.591 


1.579 


1.587 


1.603 


0.996 


1.000 


0.995 


0.900 


2.342 


1.629 


1.641 


1.622 


1.628 


1.646 


1.002 


1.009 


1.000 


0.950 


2.374 


1.699 


1.722 


1.689 


1.692 


1.710 


1.009 


1.023 


1.006 


0.980 


2.393 


1.770 


1.806 


1.756 


1.758 


1.775 


1.014 


1.037 


1.010 


Set 7 


X = 0.74 


Z = O.3Z logg 


= 14.0 


?Edd = 


1.703 keV R = 


14.80 km 


z = 0.18 


0.001 


0.303 


2.024 


2.019 


2.028 


1.351 


0.980 


0.224 


0.226 


0.223 


0.003 


0.398 


1.742 


1.731 


1.749 


1.271 


1.180 


0.426 


0.431 


0.423 


0.010 


0.538 


1.541 


1.527 


1.547 


1.269 


1.241 


0.715 


0.722 


0.711 


0.030 


0.709 


1.441 


1.444 


1.435 


1.362 


1.232 


0.910 


0.909 


0.912 


0.050 


0.805 


1.402 


1.415 


1.390 


1.386 


1.297 


0.960 


0.958 


0.964 


0.070 


0.876 


1.378 


1.397 


1.363 


1.376 


1.303 


0.976 


0.977 


0.979 


0.100 


0.957 


1.362 


1.383 


1.347 


1.365 


1.311 


0.978 


0.981 


0.980 


0.150 


1.060 


1.370 


1.384 


1.359 


1.379 


1.362 


0.975 


0.979 


0.975 


0.200 


1.139 


1.386 


1.393 


1.381 


1.407 


1.387 


0.963 


0.965 


0.962 


0.300 


1.260 


1.426 


1.424 


1.425 


1.440 


1.441 


0.978 


0.977 


0.978 


0.400 


1.354 


1.449 


1.444 


1.448 


1.464 


1.467 


0.975 


0.973 


0.975 


0.500 


1.432 


1.470 


1.464 


1.469 


1.486 


1.490 


0.977 


0.974 


0.977 


0.600 


1.498 


1.494 


1.489 


1.493 


1.507 


1.517 


0.982 


0.979 


0.981 


A *7AA 
0. /OO 


1.55 / 


1.522 


1.519 


1 COA 

1.520 


1.53 / 


1.546 


0.979 


A CiHH 
0.9/ / 


A CiHCi 

0.9/9 


0.750 


1.584 


1.546 


1.544 


1.544 


1.554 


1.570 


0.989 


0.988 


0.989 


0.800 


1.610 


1.571 


1.570 


1.568 


1.579 


1.594 


0.989 


0.988 


0.988 


0.850 


1.635 


1.605 


1.607 


1.601 


1.609 


1.627 


0.994 


0.995 


0.993 


0.900 


1.658 


1.647 


1.655 


1.642 


1.650 


1.668 


0.996 


1.001 


0.995 


0.950 


1.681 


1.720 


1.737 


1.711 


1.716 


1.737 


1.004 


1.014 


1.002 


0.980 


1.694 


1.785 


1.817 


1.772 


1.797 


1.799 


0.986 


1.004 


0.983 



Set 8 X = 0.74 Z = 0.3Z B logg = 14.3 T Edd = 2.024 keV R = 10.88 km z = 0.27 



0.001 


0.360 


1.925 


1.909 


1.933 


1.406 


1.298 


0.364 


0.371 


0.360 


0.003 


0.474 


1.685 


1.669 


1.692 


1.329 


1.269 


0.578 


0.587 


0.574 


0.010 


0.640 


1.505 


1.501 


1.504 


1.326 


1.087 


0.823 


0.825 


0.823 


0.030 


0.842 


1.386 


1.404 


1.373 


1.364 


1.281 


0.963 


0.960 


0.968 


0.050 


0.957 


1.344 


1.368 


1.329 


1.340 


1.252 


0.982 


0.982 


0.986 


0.070 


1.041 


1.330 


1.352 


1.317 


1.332 


1.304 


0.977 


0.980 


0.980 


0.100 


1.138 


1.337 


1.352 


1.328 


1.348 


1.322 


0.967 


0.970 


0.967 


0.150 


1.259 


1.364 


1.368 


1.363 


1.386 


1.364 


0.959 


0.960 


0.959 


0.200 


1.353 


1.392 


1.389 


1.393 


1.410 


1.399 


0.972 


0.970 


0.972 


0.300 


1.498 


1.420 


1.413 


1.422 


1.438 


1.434 


0.973 


0.969 


0.973 


0.400 


1.609 


1.440 


1.433 


1.441 


1.456 


1.459 


0.978 


0.973 


0.977 


0.500 


1.702 


1.460 


1.453 


1.460 


1.473 


1.484 


0.981 


0.977 


0.981 


0.600 


1.781 


1.485 


1.479 


1.484 


1.496 


1.510 


0.984 


0.980 


0.984 


0.700 


1.851 


1.520 


1.515 


1.518 


1.524 


1.545 


0.989 


0.986 


0.989 


0.750 


1.883 


1.537 


1.535 


1.535 


1.546 


1.562 


0.988 


0.986 


0.987 


0.800 


1.914 


1.564 


1.563 


1.561 


1.571 


1.587 


0.991 


0.991 


0.990 


0.850 


1.943 


1.598 


1.602 


1.594 


1.603 


1.620 


0.994 


0.996 


0.993 


0.900 


1.971 


1.645 


1.655 


1.639 


1.645 


1.664 


1.001 


1.007 


0.999 


0.950 


1.998 


1.718 


1.739 


1.709 


1.712 


1.732 


1.008 


1.020 


1.005 


0.980 


2.013 


1.776 


1.800 


1.765 


1.746 


1.786 


1.036 


1.052 


1.033 



Set 9 X = 0.74 Z = O.3Z logg = 14.6 T Edd = 2.405 keV R = 8.16 km z = 0.42 



0.001 


0.428 


1.882 


1.876 


1.883 


1.476 


1.106 


0.486 


0.489 


0.485 


0.003 


0.563 


1.651 


1.652 


1.647 


1.391 


1.206 


0.694 


0.693 


0.697 


0.010 


0.761 


1.449 


1.467 


1.436 


1.366 


1.178 


0.910 


0.900 


0.919 


0.030 


1.001 


1.317 


1.344 


1.302 


1.312 


1.209 


0.992 


0.986 


0.998 


0.050 


1.137 


1.298 


1.319 


1.288 


1.294 


1.285 


0.977 


0.976 


0.980 


0.070 


1.237 


1.310 


1.323 


1.306 


1.315 


1.301 


0.961 


0.963 


0.962 


0.100 


1.352 


1.338 


1.341 


1.340 


1.358 


1.332 


0.954 


0.954 


0.954 


0.150 


1.497 


1.372 


1.366 


1.376 


1.398 


1.370 


0.956 


0.954 


0.956 


0.200 


1.608 


1.394 


1.386 


1.397 


1.415 


1.398 


0.968 


0.964 


0.968 


0.300 


1.780 


1.415 


1.409 


1.416 


1.430 


1.430 


0.978 


0.975 


0.978 


0.400 


1.913 


1.431 


1.424 


1.430 


1.446 


1.451 


0.978 


0.975 


0.978 


0.500 


2.022 


1.456 


1.450 


1.455 


1.460 


1.481 


0.994 


0.990 


0.993 


0.600 


2.117 


1.475 


1.471 


1.474 


1.489 


1.501 


0.981 


0.978 


0.981 
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7 


r eff (keV) 


Jc,l 


fc,2 


r 

Jc,3 


hA 


fc,5 


W lfc,l 


W 2f c ,2 




0.700 


2.200 


1.514 


1.51 1 


1.513 


1.523 


1.537 


0.988 


0.986 


0.987 


0. /50 


2.238 


1.535 


1.535 


1 coo 

1.533 


1 C A 1 

1.541 


1.558 


0.992 


A AA 1 

0.991 


A AA 1 

0.991 


0.800 


2.275 


1.566 


1.567 


1.562 


1.568 


1.586 


0.996 


0.997 


0.995 


0.850 


2.309 


1.596 


1.602 


1.592 


1.598 


1.614 


0.997 


1.001 


0.996 


0.900 


2.342 


1.641 


1.654 


1.633 


1.639 


1.657 


1.002 


1.011 


1.000 


0.950 


2.374 


1.714 


1.740 


1.703 


1.707 


1.724 


1.009 


1.026 


1.006 


0.980 


2.393 


1.808 


1.851 


1.792 


1.790 


1.811 


1.021 


1.048 


1.017 


Set 10 


X = 0.74 


Z = O.1Z l0£ 


',8= 14.0 


?Edd = 


1.703 keV # = 


14.80 km 


z = 0.18 


0.001 


0.303 


1.975 


1.969 


1.979 


1.504 


1.360 


0.387 


0.390 


0.386 


0.003 


0.398 


1.767 


1.758 


1.772 


1.406 


1.188 


0.572 


0.577 


0.569 


0.010 


0.538 


1.604 


1.597 


1.606 


1.378 


1.255 


0.788 


0.791 


0.787 


0.030 


0.709 


1.495 


1.500 


1.490 


1.437 


1.261 


0.924 


0.923 


0.926 


0.050 


0.805 


1.452 


1.462 


1.445 


1.442 


1.341 


0.956 


0.956 


0.958 


0.070 


0.876 


1.430 


1.442 


1.421 


1.432 


1.345 


0.966 


0.966 


0.967 


0.100 


0.957 


1.414 


1.426 


1.406 


1.424 


1.376 


0.967 


0.969 


0.968 


0.150 


1.060 


1.412 


1.419 


1.406 


1.430 


1.400 


0.965 


0.967 


0.965 


0.200 


1.139 


1.421 


1.424 


1.418 


1.441 


1.422 


0.966 


0.967 


0.965 


0.300 


1.260 


1.440 


1.438 


1.439 


1.459 


1.453 


0.971 


0.970 


0.971 


0.400 


1.354 


1.457 


1.452 


1.456 


1.474 


1.474 


0.975 


0.972 


0.975 


0.500 


1.432 


1.474 


1.468 


1.473 


1.489 


1.493 


0.978 


0.975 


0.978 


0.600 


1.498 


1.497 


1.492 


1.496 


1.510 


1.519 


0.982 


0.979 


0.981 


A -7AA 
0. /OO 


1.55 / 


1.522 


1.519 


1 COA 

1.520 


1.539 


1.546 


0.977 


A ATC 
0.9/5 


0.9/ / 


0.750 


1.584 


1.551 


1.548 


1.550 


1.558 


1.574 


0.990 


0.988 


0.990 


0.800 


1.610 


1.574 


1.574 


1.571 


1.580 


1.597 


0.992 


0.992 


0.991 


0.850 


1.635 


1.608 


1.610 


1.604 


1.613 


1.629 


0.992 


0.993 


0.991 


0.900 


1.658 


1.651 


1.659 


1.646 


1.655 


1.672 


0.995 


0.999 


0.994 


0.950 


1.681 


1.730 


1.748 


1.721 


1.724 


1.746 


1.006 


1.017 


1.004 


0.980 


1.694 


1.787 


1.811 


1.778 


1.763 


1.802 


1.031 


1.044 


1.028 



Set 11 X = 0.74 Z = O.1Z log g = 14.3 T Edd = 2.024 keV R = 10.88 km z = 0.27 



0.001 


0.360 


1.910 


1.900 


1.915 


1.537 


1.302 


0.519 


0.525 


0.516 


0.003 


0.474 


1.719 


1.711 


1.722 


1.445 


1.283 


0.693 


0.698 


0.691 


0.010 


0.640 


1.560 


1.560 


1.556 


1.416 


1.275 


0.866 


0.866 


0.868 


0.030 


0.842 


1.444 


1.456 


1.435 


1.424 


1.314 


0.956 


0.955 


0.959 


0.050 


0.957 


1.407 


1.420 


1.398 


1.404 


1.319 


0.967 


0.968 


0.969 


0.070 


1.041 


1.393 


1.404 


1.386 


1.400 


1.354 


0.966 


0.967 


0.966 


0.100 


1.138 


1.390 


1.397 


1.387 


1.407 


1.369 


0.962 


0.964 


0.962 


0.150 


1.259 


1.402 


1.402 


1.401 


1.424 


1.399 


0.962 


0.962 


0.962 


0.200 


1.353 


1.414 


1.410 


1.414 


1.435 


1.419 


0.967 


0.965 


0.967 


0.300 


1.498 


1.430 


1.424 


1.430 


1.448 


1.444 


0.973 


0.970 


0.973 


0.400 


1.609 


1.444 


1.437 


1.444 


1.460 


1.463 


0.977 


0.973 


0.977 


0.500 


1.702 


1.463 


1.455 


1.462 


1.476 


1.485 


0.981 


0.976 


0.981 


0.600 


1.781 


1.488 


1.481 


1.487 


1.499 


1.512 


0.984 


0.980 


0.984 


0.700 


1.851 


1.522 


1.517 


1.521 


1.532 


1.546 


0.987 


0.983 


0.987 


0.750 


1.883 


1.540 


1.537 


1.539 


1.553 


1.563 


0.983 


0.981 


0.983 


0.800 


1.914 


1.569 


1.569 


1.566 


1.574 


1.592 


0.994 


0.994 


0.993 


0.850 


1.943 


1.603 


1.607 


1.599 


1.607 


1.625 


0.995 


0.997 


0.994 


0.900 


1.971 


1.649 


1.659 


1.643 


1.650 


1.667 


0.999 


1.005 


0.997 


0.950 


1.998 


1.728 


1.750 


1.718 


1.724 


1.741 


1.006 


1.019 


1.003 


0.980 


2.013 


1.789 


1.816 


1.778 


1.758 


1.799 


1.038 


1.055 


1.035 


Set 12 


X = 0.74 


Z = O.1Z logg=14.6 


?Edd 


= 2.405 keV R 


= 8.16 km 


z = 0.42 


0.001 


0.428 


1.868 


1.866 


1.868 


1.581 


1.377 


0.623 


0.625 


0.623 


0.003 


0.563 


1.675 


1.679 


1.673 


1.485 


1.217 


0.784 


0.782 


0.786 


0.010 


0.761 


1.503 


1.517 


1.495 


1.431 


1.193 


0.925 


0.919 


0.930 


0.030 


1.001 


1.394 


1.408 


1.386 


1.381 


1.330 


0.969 


0.968 


0.972 


0.050 


1.137 


1.373 


1.383 


1.369 


1.374 


1.331 


0.964 


0.964 


0.964 


0.070 


1.237 


1.374 


1.379 


1.374 


1.386 


1.351 


0.959 


0.960 


0.959 


0.100 


1.352 


1.385 


1.384 


1.386 


1.405 


1.373 


0.959 


0.959 


0.959 


0.150 


1.497 


1.399 


1.395 


1.400 


1.420 


1.399 


0.965 


0.963 


0.965 


0.200 


1.608 


1.408 


1.403 


1.409 


1.427 


1.414 


0.972 


0.970 


0.972 


0.300 


1.780 


1.419 


1.414 


1.419 


1.434 


1.435 


0.977 


0.975 


0.977 


0.400 


1.913 


1.433 


1.427 


1.433 


1.448 


1.453 


0.980 


0.976 


0.979 
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7 


r eff (keV) 


Jc,l 


Jc,2 


/c,3 


JcA 


fc,5 


W lfc,l 


W 2fc,2 




0.500 




2.022 


1.454 


1.448 


1.453 


1.466 


I All 


0.983 


0.979 


0.983 


0.600 




2.117 


1.478 


1.473 


1.477 


1.491 


1.502 


0.982 


0.979 


0.982 


A 1AA 

0.700 




2.200 


1.516 


1 C 1 1 

1.513 


1.515 


1.524 


1.539 


A AO O 

0.988 


A AO/C 

0.986 


A AO O 

0.988 


0. /50 




2.238 


1.537 


1.537 


1.534 


1 C A 1 

1.543 


1.559 


0.992 


A AA 1 

0.991 


A AA 1 

0.991 


0.800 




2.275 


1.567 


1.569 


1.564 


1.569 


1.588 


0.997 


0.998 


0.996 


0.850 




2.309 


1.600 


1.605 


1.595 


1.602 


1.617 


0.997 


1.000 


0.996 


0.900 




2.342 


1.646 


1.659 


1.639 


1.647 


1.662 


0.999 


1.007 


0.997 


0.950 




2.374 


1.718 


1.745 


1.707 


1.712 


1.728 


1.008 


1.025 


1.005 


0.980 




2.393 


1.795 


1.839 


1.778 


1.798 


1.798 


0.997 


1.024 


0.992 


Set 13 


X 


= 0.74 


Z = 0.01Z o lo g<? =14.0 


?Edd = 


1.703 keV R 


= 14.80 km 


z = 0.18 


0.001 




0.303 


1.949 


1.945 


1.951 


1.713 


1.521 


0.668 


0.672 


0.666 


0.003 




0.398 


1.780 


I. Ill 


1.782 


1.591 


1.427 


0.780 


0.783 


0.779 


0.010 




0.538 


1.644 


1.642 


1.644 


1.512 


1.305 


0.880 


0.881 


0.880 


0.030 




0.709 


1.550 


1.551 


1.548 


1.512 


1.361 


0.938 


0.937 


0.938 


0.050 




0.805 


1.515 


1.517 


1.512 


1.512 


1.398 


0.953 


0.953 


0.954 


0.070 




0.876 


1.495 


1.497 


1.492 


1.503 


1.412 


0.960 


0.960 


0.960 


0.100 




0.957 


1.477 


1.480 


1.474 


1.492 


1.438 


0.965 


0.966 


0.965 


0.150 




1.060 


1.462 


1.463 


1.459 


1.480 


1.449 


0.969 


0.969 


0.969 


0.200 




1.139 


1.456 


1.456 


1.453 


1.474 


1.456 


0.971 


0.971 


0.971 


0.300 




1.260 


1.455 


1.452 


1.453 


1.472 


1.465 


0.974 


0.973 


0.974 


0.400 




1.354 


1.462 


1.458 


1.460 


1.478 


1.479 


0.977 


0.974 


0.976 


0.500 




1.432 


1.476 


1.471 


1.475 


1.490 


1.495 


0.979 


0.976 


0.978 


0.600 




1.498 


1.497 


1.492 


1.496 


1.511 


1.520 


0.981 


0.978 


0.981 


A TAA 

0.700 




1.557 


1.529 


1.524 


1.528 


1.542 


1.552 


0.982 


A ATA 

0.979 


A flOO 

0.982 


0. /50 




1.584 


1.549 


1.547 


1.547 


1 <^A 

1.560 


1.572 


0.986 


A AO A 

0.984 


A AO C 

0.985 


0.800 




1.610 


1.575 


1.574 


1.573 


1.586 


1.597 


0.986 


0.985 


0.986 


0.850 




1.635 


1.608 


1.611 


1.604 


1.613 


1.630 


0.993 


0.994 


0.992 


0.900 




1.658 


1.654 


1.661 


1.649 


1.659 


1.674 


0.994 


0.998 


0.993 


0.950 




1.681 


1.732 


1.750 


1.724 


1.730 


1.748 


1.004 


1.014 


1.002 


0.980 




1.694 


1.803 


1.836 


1.789 


1.815 


1.814 


0.986 


1.006 


0.983 


Set 14 


X 


= 0.74 


Z = 0.01Z o log 


8 = 14.3 


?Edd 


2.024 keV R 


= 10.88 km 


z = 0.27 


0.001 




0.360 


1.900 


1.896 


1.902 


1 .700 


1.520 


0.733 


0.737 


0.731 


0.003 




0.474 


1.739 


1.735 


1.740 


1.580 


1.423 


0.830 


0.832 


0.829 


0.010 




0.640 


1.606 


1.606 


1.606 


1.509 


1.348 


0.909 


0.909 


0.909 


0.030 




0.842 


1.516 


1.518 


1.513 


1.502 


1.378 


0.951 


0.951 


0.951 


0.050 




0.957 


1.483 


1.485 


1.481 


1.491 


1.406 


0.961 


0.961 


0.961 


0.070 




1.041 


1.466 


1.468 


1.463 


1.479 


1.418 


0.965 


0.966 


0.965 


0.100 




1.138 


1.451 


1.452 


1.449 


1.468 


1.429 


0.968 


0.968 


0.968 


0.150 




1.259 


1.440 


1.439 


1.438 


1.457 


1.436 


0.971 


0.971 


0.971 


0.200 




1.353 


1.435 


1.433 


1.434 


1.453 


1.441 


0.974 


0.973 


0.973 


0.300 




1.498 


1.437 


1.432 


1.436 


1.452 


1.451 


0.977 


0.974 


0.976 


0.400 




1.609 


1.446 


1.440 


1.446 


1.462 


1.465 


0.978 


0.974 


0.978 


0.500 




1.702 


1.459 


1.452 


1.458 


1.479 


1.482 


0.971 


0.967 


0.971 


0.600 




1.781 


1.494 


1.486 


1.493 


1.498 


1.517 


0.994 


0.990 


0.994 


0.700 




1.851 


1.525 


1.518 


1.524 


1.534 


1.548 


0.988 


0.983 


0.987 


0. /50 




1.883 


1.546 


1.542 


1.545 


1.555 


1.568 


0.988 


A AO C 

0.985 


A AO O 

0.988 


0.800 




1.914 


1.571 


1.570 


1.569 


1.579 


1.593 


0.990 


0.989 


0.990 


0.850 




1.943 


1.604 


1.607 


1.600 


1.609 


1.624 


0.993 


0.995 


0.992 


0.900 




1.971 


1.651 


1.661 


1.645 


1.652 


1.669 


0.999 


1.005 


0.997 


0.950 




1.998 


1.728 


1.749 


1.718 


1.721 


1.741 


1.008 


1.021 


1.006 


0.980 




2.013 


1.775 


1.800 


1.765 


1.746 


1.786 


1.037 


1.052 


1.034 


Set 15 


X 


= 0.74 


Z = 0.01Z o loj 


lg = 14.6 r Edd = 


: 2.405 keV R 


= 8.16 km 


z = 0.42 


0.001 




0.428 


1.861 


1.858 


1.861 


1.695 


1.427 


0.778 


0.780 


0.777 


0.003 




0.563 


1.702 


1.702 


1.701 


1.576 


1.427 


0.864 


0.864 


0.864 


A A 1 A 
0.010 




0.761 


1.571 


1.573 


1.569 


1.497 


1.378 


0.929 


A AO O 

0.928 


A AO A 

0.930 


0.030 




1.001 


1.485 


1.488 


1.482 


1.476 


1.387 


0.959 


0.959 


0.959 


0.050 




1.137 


1.456 


1.459 


1.454 


1.465 


1.402 


0.965 


0.965 


0.965 


0.070 




1.237 


1.442 


1.443 


1.440 


1.455 


1.412 


0.968 


0.968 


0.968 


0.100 




1.352 


1.430 


1.430 


1.428 


1.445 


1.418 


0.971 


0.971 


0.971 


0.150 




1.497 


1.421 


1.420 


1.419 


1.437 


1.422 


0.974 


0.974 


0.974 


0.200 




1.608 


1.418 


1.416 


1.416 


1.433 


1.427 


0.976 


0.975 


0.976 
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7 


7e ff (keV) 


£ 

Jc,l 


£ 

fc,2 


£ 

Jc,3 


£ 

JcA 


£ 

fc,5 


W lfc,l 


W 2f c ,2 


W lfc.3 


0.300 


1.780 


1.421 


1.417 


1.420 


1.436 


1.438 


0.978 


0.976 


0.978 


0.400 


1.913 


1.434 


1.429 


1.434 


1.448 


1.454 


0.980 


0.976 


0.980 


0.500 


2.022 


1.455 


1.448 


1.455 


1.467 


1.477 


0.983 


0.979 


0.983 


0.600 


2.117 


1.484 


1.478 


1.484 


1.493 


1.507 


0.988 


0.983 


0.988 


0.700 


2.200 


1.518 


1.514 


1.516 


1.526 


1.540 


0.989 


0.987 


0.989 


0.750 


2.238 


1 C A A 

1.540 


1 CIO 

1.538 


1.538 


1.548 


1.562 


0.989 


0.988 


a no o 
0.988 


0.800 


2.275 


1.570 


1.571 


1.567 


1.573 


1.590 


0.996 


0.996 


0.995 


0.850 


2.309 


1.601 


1.607 


1.596 


1.603 


1.619 


0.998 


1.002 


0.997 


0.900 


2.342 


1.649 


1.663 


1.642 


1.647 


1.665 


1.003 


1.011 


1.001 


0.950 


2.374 


1.724 


1.751 


1.713 


1.717 


1.734 


1.009 


1.025 


1.006 


0.980 


2.393 


1.785 


1.816 


1.773 


1.750 


1.791 


1.042 


1.062 


1.039 


Set 16 7 = 1 


Z = 


logg = 


14.0 r Edd = 1.955 keV 


R = 14.80 km 


z = 0.18 


0.001 


0.348 


1.967 


1.956 


1.973 


1.781 


1.609 


0.793 


0.801 


0.788 


0.003 


0.457 


1.785 


1.774 


1.791 


1.637 


1.517 


0.862 


0.869 


0.858 


0.010 


0.618 


1.641 


1.628 


1.646 


1.552 


1.423 


0.914 


0.917 


0.912 


0.030 


0.814 


1.545 


1.534 


1.550 


1.550 


1.375 


0.946 


0.945 


0.945 


0.050 


0.924 


1.509 


1.498 


1.513 


1.529 


1.379 


0.957 


0.954 


0.957 


0.070 


1.005 


1.488 


1.477 


1.491 


1.508 


1.390 


0.962 


0.959 


0.962 


0.100 


1.099 


1.467 


1.457 


1.470 


1.487 


1.405 


0.967 


0.964 


0.968 


0.150 


1.217 


1.446 


1.435 


1.448 


1.464 


1.419 


0.972 


0.967 


0.972 


0.200 


1.307 


1.433 


1.421 


1.435 


1.450 


1.424 


0.974 


0.968 


0.974 


0.300 


1.447 


1.419 


1.405 


1.421 


1.433 


1.428 


0.976 


0.968 


0.977 


0.400 


1.555 


1.415 


1.400 


1.417 


1.430 


1.432 


0.977 


0.968 


0.977 


0.500 


1.644 


1.419 


1.402 


1.421 


1.434 


1.441 


0.978 


0.967 


0.978 


0.600 


1.720 


1.430 


1.412 


1.433 


1.445 


1.455 


0.979 


0.968 


0.980 


0.700 


1.788 


1 A C 1 

1.451 


1 All 

1.433 


1.454 


1.465 


1.478 


0.980 


0.968 


a no 1 
0.981 


0.750 


1.819 


1.468 


1.451 


1.470 


1.480 


1.494 


0.983 


0.971 


0.983 


0.800 


1.849 


1.487 


1.472 


1.489 


1.500 


1.514 


0.982 


0.972 


0.983 


0.850 


1.877 


1.518 


1.504 


1.520 


1.529 


1.543 


0.985 


0.976 


0.985 


0.900 


1.904 


1.556 


1.548 


1.556 


1.569 


1.579 


0.984 


0.978 


0.983 


0.950 


1.930 


1.629 


1.634 


1.624 


1.632 


1.649 


0.996 


0.999 


0.995 


0.980 


1.945 


1.686 


1.688 


1.682 


1.659 


1.703 


1.033 


1.034 


1.032 



Set 17 7=1 Z = logg = 14.3 T Edd = 2.323 keV R = 10.88 km z = 0.27 



0.001 


0.413 


1.924 


1.913 


1.929 


1.751 


1.598 


0.821 


0.828 


0.817 


0.003 


0.544 


1.748 


1.736 


1.754 


1.617 


1.508 


0.884 


0.890 


0.881 


0.010 


0.735 


1.609 


1.597 


1.614 


1.555 


1.418 


0.929 


0.931 


0.928 


0.030 


0.967 


1.517 


1.505 


1.521 


1.530 


1.378 


0.956 


0.954 


0.955 


0.050 


1.099 


1.482 


1.470 


1.485 


1.503 


1.384 


0.964 


0.961 


0.964 


0.070 


1.195 


1.461 


1.450 


1.463 


1.480 


1.393 


0.969 


0.965 


0.969 


0.100 


1.306 


1.440 


1.429 


1.443 


1.458 


1.403 


0.973 


0.968 


0.973 


0.150 


1.446 


1.420 


1.408 


1.422 


1.436 


1.409 


0.976 


0.970 


0.976 


0.200 


1.554 


1.408 


1.394 


1.410 


1.423 


1.411 


0.977 


0.970 


0.977 


0.300 


1.719 


1.396 


1.380 


1.399 


1.411 


1.412 


0.978 


0.968 


0.979 


0.400 


1.848 


1.395 


1.377 


1.399 


1.410 


1.417 


0.979 


0.967 


0.979 


0.500 


1.954 


1.403 


1.382 


1.406 


1.416 


1.428 


0.980 


0.966 


0.980 


0.600 


2.045 


1.417 


1.396 


1.421 


1.431 


1.444 


0.980 


0.966 


0.981 


0.700 


2.125 


1.443 


1.423 


1.446 


1.452 


1.471 


0.987 


0.972 


0.988 


0.750 


2.162 


1.461 


1.442 


1.464 


1.466 


1.489 


0.993 


0.979 


0.993 


0.800 


2.197 


1.480 


1.462 


1.483 


1.497 


1.506 


0.978 


0.965 


0.979 


0.850 


2.231 


1.514 


1.500 


1.515 


1.521 


1.539 


0.989 


0.979 


0.989 


0.900 


2.263 


1.559 


1.551 


1.558 


1.562 


1.580 


0.995 


0.990 


0.995 


0.950 


2.294 


1.630 


1.636 


1.625 


1.630 


1.646 


1.001 


1.005 


0.999 


0.980 


2.312 


1.690 


1.699 


1.683 


1.664 


1.702 


1.032 


1.039 


1.030 


Set 18 7=1 


Z = 


logg = 


14.6 


r Edd = 2.761 keV 


R = 8 


.16 km 


z = 0.42 


0.001 


0.491 


1.893 


1.884 


1.897 


1.734 


1.590 


0.838 


0.844 


0.834 


0.003 


0.646 


1.721 


1.712 


1.726 


1.605 


1.501 


0.897 


0.902 


0.895 


0.010 


0.873 


1.584 


1.575 


1.588 


1.546 


1.412 


0.939 


0.940 


0.938 


0.030 


1.149 


1.492 


1.484 


1.495 


1.506 


1.380 


0.963 


0.961 


0.963 


0.050 


1.306 


1.457 


1.449 


1.459 


1.473 


1.386 


0.970 


0.967 


0.970 


0.070 


1.420 


1.437 


1.429 


1.438 


1.452 


1.392 


0.973 


0.971 


0.974 


0.100 


1.553 


1.417 


1.408 


1.418 


1.432 


1.396 


0.976 


0.973 


0.976 
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7 


r eff (keV) 


r 

Jc,l 


fc,2 


r 

Jc,3 


JcA 


fc,5 


W lfc,l 


W 2fc,2 


W lfc.3 


0.150 


1.718 


1.397 


1.388 


1.399 


1.411 


1.396 


0.978 


0.973 


0.979 


0.200 


1.847 


1.387 


1.375 


1.388 


1.400 


1.394 


0.979 


0.973 


0.979 


0.300 


2.043 


1.378 


1.363 


1.380 


1.392 


1.396 


0.980 


0.970 


0.980 


0.400 


2.196 


1.380 


1.363 


1.384 


1.394 


1.402 


0.980 


0.968 


0.981 


0.500 


2.322 


1.390 


1.370 


1.394 


1.404 


1.415 


0.979 


0.965 


0.980 


0.600 


2.430 


1.409 


1.389 


1.412 


1.421 


1.435 


0.982 


0.968 


0.983 


0. /OO 


2.526 


1.43 / 


1 A 1 *7 

1 .4 1 / 


1 A A A 

1.440 


1 A <A 

1.450 


1 A 

1.462 


A AO 1 


0.966 


A AOO 

0.9o2 


0.750 


2.570 


1.457 


1.439 


1.460 


1.466 


1.483 


0.987 


0.974 


0.988 


0.800 


2.611 


1.481 


1.466 


1.482 


1.488 


1.505 


0.990 


0.979 


0.991 


0.850 


2.651 


1.513 


1.502 


1.514 


1.519 


1.535 


0.992 


0.984 


0.992 


0.900 


2.689 


1.557 


1.552 


1.555 


1.560 


1.576 


0.995 


0.992 


0.994 


0.950 


2.726 


1.629 


1.640 


1.622 


1.625 


1.641 


1.005 


1.012 


1.002 


0.980 


2.747 


1.678 


1.688 


1.671 


1.646 


1.688 


1.039 


1.047 


1.037 



